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^ '. Abstract 

The "positive square" of any tensor is presented in a universal and unified man- 
■ ner, valid in Lorentzian manifolds of arbitrary dimension, and independently of any 

(anti)-symmetry properties of the tensor. For rank-m tensors, the positive square 
has rank 2m. Positive here means future, that is to say, satisfying the dominant 
property. The standard energy-momentum and super-energy tensors are recovered 
as appropriate parts of the general square. A richer structure of principal null 
directions arises. 
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PACS: 04.20.-q, 04.20.Cv, 02.40.-k 



The energy-momentum tensors of a 2-form (such as the electromagnetic field) or of 
a 1-form (the gradient of the scalar field), or the Bel [3] and Bel- Robinson [3| tensors 
('super-energy tensor' of the gravitational field, see also JO]), or the super-energy tensors 
of the electromagnetic jTH ES| and scalar fields [2H1 ESI Ell- etc., have either different 
expressions or different number of indices, or both. A unified, algebraic, perspective of 
all these tensors was presented in [23] for Lorentzian manifolds of any dimension, where 
a general method to construct the positive square 1 of any tensor, called generically its 
'superenergy' tensor, was put forward. This method recovered the standard expressions for 
the energy-momentum and super-energy tensors of the physical fields, keeping therefore 
its diversity in rank and structure, and allows us to understand that these differences arise 
because of, and depend on, (i) the number of anti-symmetric sets into which the indices 
of the fields fall, and (ii) the number of indices within each of these antisymmetric blocks. 

In this short Note, I would like to call attention to a very useful positive square of a 
tensor that does not seem to have been noticed before. It is the most general — under 
some conditions — future tensor quadratic on the given tensor. 

1 By "positive square" of a tensor I mean a future tensor quadratic in the original one. Future tensors 
are those satisfying the dominant property, see (JSJ below; see section 4 in section 2 in [7], section 
2.1 in JT§|, section 7 in 23 , and references therein for further details. 
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Let (V, g) be any Lorentzian manifold of arbitrary dimension n, with metric tensor g 
and signature (— , +,...,+). The basic structure to be used is given by the very simple 
tensor 

G^* = + g^sf* (1) 

which has the properties 

G\V Q/3 = ^(Ajtt) = G^xn, G p pa/ 3 = G aj 3 P p = (2 — n) g a p, G p apl 3 = ng aj3 , (2) 

and when considered as an endomorphism on rank-2 tensors acts as follows: 

G\fj, al3 t a p = 2t(X/i) — t p P g\n ■ 

The 'universal (super)-energy operator' is simply given by the m th power of G, that is to 
say, the tensor product ® m G, where the number of multiplications m is the rank of the 
tensor to be squared. With indices this operator reads 




Then, given any tensor A fXl ,__ olm one can define its future square by letting this operator 
act on the tensor product A® A 

£\ini...\ m n m {A} = - \ Y[G(x if i,) ai % J A ai ... am Aj3 l ...p m = 

= 2^Ai/ tl Ql/31 • • ■ G\ m p m am P m A ai _ am Ap 1 _p m . (3) 

Observe that this definition is independent of the dimension n. In order to distinguish this 
generalization from its predecessors — the energy-momentum and super-energy tensors — 
I propose to call the tensor (J3J) the "mathematical energy tensor" of A. 2 Note that one 
may let the operator act on two different tensors A^ ... Mm and ^ as long as they have 
the same rank m; an equivalent way of defining this operation is (suppressing indices) 

S{A, B} = 28 ~ \S{A} - U{B} 

so that £{A, B} = £{B, A} and £{A, A} = £{A}. 

Now, let me enumerate the more important properties of the tensors (J3J). 

(i) The mathematical energy tensor of a rank-m tensor has rank 2m. From (J2J) and 
(HH) one derives 

^Ai/ii...A m /i m {A B} = £(\ lf M 1 )...(\ mf i m ){A, B} , 

2 Other possible names are the "causal squares" of A for ±£{A}, and in particular the "future square" 
of A for the plus sign. 
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and similarly for other traces on any pair A^. 

(ii) As is obvious from (j3J), if A ai _ am is symmetric, or antisymmetric, with respect to 
two indices a« and oej, then £\ llJ , 1 ...\ ml j m {A} is symmetric on the interchange of the corre- 
sponding pairs XiHi and Xjfij. Similar properties can be listed for £Ai^i...A m /i m {A B} if the 
two tensors have the (anti)-symmetries on corresponding indices. Thus, if A ai „. ai _ a ._ am = 
eiA*i...cy...a<...a m and /£,....„....<>,.....,.. = ^2B ai ...a j ...a i ...a m were e\ = e\ = 1, then 

{A,B} 

I'll 

{A,B}. 

(iii) Observe that (jSJ) coincides with the basic superenergy tensor of any tensor A ai _ am 
with no antisymmetries. To see this, let me recall that the superenergy construction is 
based on the following fact [21] |23 E] : there is a precise and unique way to consider 
any tensor A fJ/1 ___ fJim as an r-fold form, that is to say, as an element of A™ 1 eg) . . . eg A nr 
where A p is the set of p-forms. The well defined natural number r is called the form- 
structure number of A^...^ [13]; obviously, r < m. Moreover, the set of r natural 
numbers n\, . . . ,n r is univocally defined, and each ny is called the T-th block rank. It is 
trivial that J2r=i n x = m- Tensors seen in this way are called r-fold (rii, . . . , n r )-forms 
and denoted by ^4[ ni ]...[ nr ] [22] • 

Given any tensor A^^ ^ with form structure number r, its basic superenergy tensor 
TAi^i...A r ^t r {^4} has 2r indices, it is symmetric on each of the r pairs, and can be defined 
by |23 EI]: 3 

T\ \ $A\ = -E\ \ a 1 ...a ni ...T 1 ...T nr 'X j{p 1 ...p ni ...v 1 ...Vn r (A) 

-<- A 1 fl 1 ...\ r fl r l^J 2 ■^ll J -l---^rlJ.r p 1 ...p ni ...V\...V Hr ■ ■ ■°'n 1 ■ ■ -TL . . ,T„ r V V 

where 

P. . cr 1 ...a ni ...T 1 ...T nr — _ S(T 2 -cr ni (n^i n La "!^, \ -x ■ ■ ■ 

± ~'Aip,i...\rllrp 1 ... Pni ...ui...v nr — _ -Qj u p2...pn 1 ^ (Ai Jl\ P1 ^ J 

1 „ - 1 



and $vl'"jf£ = P'^fti ' ' ' * s ^ ne Kronecker symbol of order p. 

It follows that, for the case when A /Jil ___ IXm has no antisymmetries, that is to say, when 
r = m, the two tensors (|3J) and coincide: 

If r = m => S Xltll ...\ m ^ m {A} = T Al/il ... AmMm {A} . 

Thus, for those familiar with superenergy tensors, the rule in order to compute (jSJ) 
for general tensors is very simple: just ignore the antisymmetries of the tensor A, its 
form structure number and block ranks, and build its super-energy tensor as an m-fold 
(1,1,. . . ,l)-form. 

3 Here, A^...^ is the permuted version of ^4 Ml ... Mm so that its first m indices are those precisely in 
[m], the next ri2 indices are those in [712], and so on. 
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(iv) The most important property is that the mathematical energy tensor Q is a 
future tensor, see [3 UHl d] . This means that it always satisfies the dominant property 
for any collection of future-directed vectors {ui, Vi, . . . , u m , v m } 

Sx lfll ...x mfim {A}u^ ■ --u x -v^ > . (5) 

This property permits the finding of estimates on the growth of the tensor field A (if it 
is subject to some differential field equations) and thereby has been essential in many 
important applications of the energy-momentum and super-energy tensors, such as the 
causal propagation of fields [201 El E] , the stability of spacetimes [T31 [21] , getting hyper- 
bolizations of many physical fields including gravity [13 EH d El U21 121] , the existence 
of global solutions to the gravitational Cauchy problem [2TJ E] , and the propagation of 
discontinuities [221 125] . 

The basic super-energy tensor (@J) of A is the most general rank-2r tensor — modulo 
index permutations, see section 5 in [23] — , symmetric by pairs and quadratic on A satis- 
fying the dominant property As a matter of fact, the tensor (J3J) is analogously the most 
general tensor, modulo index permutations, quadratic on A (and without extra g's, that 
is, with rank 2m) satisfying the dominant property — and symmetric on each Aj//;-pair 
(see footnote El) 

The dominant property can be characterized by several different but equivalent con- 
ditions, see [23 HH] . Among them, one can cite (a) that the inequality (0) is strict if 
all the future-directed vectors {ux,Vi, . . . ,u m ,v m } are timelike (if A is non-vanishing, of 
course); (b) that the vectors 

' — C Mi-AmMml^J^l ■■■ U m V m 

are future directed for any arbitrary set of future-directed vectors {vi, . . . ,u m ,v m }, and 
analogously by leaving any other single index free; or (iii) that the totally timelike com- 
ponent 'dominates' the tensor, meaning that in any orthonormal basis {e*o, . . . , e n _i} with 
eo timelike future-directed, 

£oO...Oo{^4} ^ |£XiA«i...AmMm {^4}l (6) 

for all possible values of the indices. 

To prove that the dominant property holds for (JSJ) one can just use any of the available 
proofs 4 for the superenergy tensors [23I2S], apply it to an m-fold (1,. . . ,l)-form, and note 
that any possible antisymmetries of the tensor A do not invalidate the proofs. As the 
property holds for generic tensors A, then nothing changes if the tensor A happens to 
have some antisymmetries. 

4 Here is where the symmetrization on each Ai/z^-pair is needed, as there is no available proof of the 
dominant property in general dimension without this restriction. However, in four dimensions, this 
can be generalized by using spinors and the reasonings in Therefore, it seems plausible that these 
symmetrizations can be avoided in general dimension. Observe that the superenergy density of the next 
point (v) is not affected by this at all. 
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(v) The previous point concerning the validity of the general proof of the dominant 
property is actually related to the fact that the generic 'electric-magnetic' decomposi- 
tion of any tensor ignoring its antisymmetries is essentially the same as the canonical 
'electric-magnetic' decomposition (see [21 "2H] and section 2 in [2*5] ) in which the form 
structure number and block ranks play an important role. The reason for this is that 
the canonical 'electric-magnetic' decomposition is nothing but the generic one with some 
'electric-magnetic' parts — those corresponding to more than one 'electric', or more than 
one 'magnetic', part within one antisymmetric block — vanishing. 

Using this fact, one can easily derive the fact that the totally timelike component 
arising in (J""} is, up to a multiplicative factor, the so-called superenergy density of the 
tensor A relative to the unit timelike vector eo. This is denoted by Wa(co), it is defined 
by [22] 

W A (e ) = T Aim ... Ar ^{A}e Al e^ . . . 

and it is equal to half the sum of the positive squares of all the canonical electric-magnetic 
parts of the tensor A m ...^ m [25 . The supernergy density is also equal to half the sum of 
the squares of all the independent components of ^4 Ml ... Mm in any orthonormal basis {e M } 
containing e$ (this is the typical mathematical 'energy' of the tensor A^ ... Atm ): 5 



\ ( r \ \ 71—1 



W A (e )=T 00 ... {A} = -[U—) £ IA 



|2 



T=l nr -J m,...,lMn=0 



It follows then easily that 



n-1 



£oo...oo{A} = \ £ lA^..^ 2 = ( II n r\ ) W A (e ) . 

/ ^i,...,/i m =0 \T=1 / 



Therefore, from this point of view, the new tensor (["")) and the superenergy tensor (["")) are 
essentially equivalent. 

(vi) Point (v) leads to the idea that the superenergy tensor of A is in fact contained 
as an essential part of the general tensor ([""J), and some symmetric part of the latter is 
probably nothing but the former supplemented with the necessary g's. This is true, and 
the proof is as follows. Observe that the first summands in the explicit expression of 
£\ lfll ...\ m n m {A} are always of the type A\ 1 ^\ m A fMl fMm and the symmetrizations of this on 
each Aj/ij pair. There are 2 m_1 terms of this type. The next summands are of type 



9\ 1 ii 1 Ap\ 2 ...\ m A p 



and the symmetrizations of this on each Xi/j,i pair, together with the similar ones with an 
explicit g\ tfll and the corresponding contraction on the i th index of A. There are a total of 
m2 m ~ 2 summands of this type. And so on by extracting explicit metrics and contracting 
on the corresponding indices of A <g> A. 

Now, suppose that A\ ± x m has some antisymmetries. Then, the symmetric part of the 
first type of terms is identically zero, and therefore its contribution to the mathematical 

5 There is a typo in Property 3.5 of "2S|: a missing factor I1t=i nf- 
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energy vanishes. This is also the case for the terms of the second type corresponding to a 
g Xitii if the indices of A Xl _ Xm with \ excluded still have antisymmetries. And so on for the 
rest of the summands until one arrives at terms with a number of explicit metrics given 
by Z)x=i( n T — 1) where r is the form structure number of A. Therefore, the following 
formula hods in general 

^ lft ...a„ l/3ni ) (7l5l - 7 " 2<Sn2) (A lW ...A„ rMnr ){i} = [fii-ir- 1 ny\\ X 

9(ofh ■ ■■g ani p n T ai p l) ^ s \.. {Xl , 1 {A}g^ ■ ..g^u 2 ) . .. gx2fi2 . .. gXnrfinr) . ( 7 ) 

Thus, the standard basic superenergy tensor is fully contained in (jHJ), and can be obtained 
from © by means of the important formula (J7|). Observe, nevertheless, that the mathe- 
matical energy tensor Q has more information than as the positive property (0) for 
different future-pointing vectors implies. 

(vii) The principal directions of a future tensor are defined by the causal vectors k 
such that the contraction of the tensor on all its indices with k vanishes, see Definition 
A. 2 in [THj. As remarked above, the inequality (0) is strict if all the vectors are timelike. 
Therefore, due to the general properties of causal tensors [3 [TH] , the principal directions 
are necessarily null. The principal directions of the tensor Q are thus given in general 
by the null vectors k such that 

k^k^...k x -k^£ Xl , 1 ... Xm , m {A} = 0. (8) 

Of course, if the tensor A has no antisymmetries, i.e., when r = m, then these are simply 
the principal directions of the superenergy tensor of A, which are in this case the null 
directions k such that k Xl ■ ■ ■ k m A Xl ...x m = 0. However, if r < m, that is to say, if the 
tensor A has at least two antisymmetric indices, then every null vector is a principal 
direction of £{A}: 

lir<m => k x ^^---k x ™k^S Xllll ... XmlIrn {A} = Vfc = 0). 

Fortunately, this does not mean that the information on the principal null directions 
of A (see e.g. Definition 2 in [23J), which coincide with the principal directions of its 
superenergy tensor and are thus characterized by 

k^k^---k x ^k^T Xl , 1 ... Xr ^{A} = 0, 

has been lost. This can also be deduced from point (vi). To recover this information, one 
only has to start suppressing copies of k from formula (jSjl. As a matter of fact, if one 
does this in order and systematically, the information about the form structure number 
and block ranks of A can be retrieved. 

The procedure is based on the fundamental relation ((7|) and is as follows. Start by 
suppressing one copy of k, and do this for each Aj/ij-pah- That is to say, for each i = 
1, . . . , m, compute 

fcAifcw ■ ■ ■ k x < ■ ■ ■ k^k^£ Xltll ... Xm ... Xmlhn {A} 
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(here k^ is absent.) If for a given i there is at least one null k such that the previous 
expression is non-zero, then there exists exactly one j ^ i such that the same happens 
by removing only the k^ . In this case the only antisymmetric indices in A fJil ___ IXm are the 
couple fa and fa:. This is the end of the story in this case, r = m — 1 and there are 
m — 1 blocks with one index and one block with rank 2. If on the contrary, the previous 
displayed expression vanishes for all null k and all values of % = l,...,m, then either 
r < in — 1 or r = m — 1 and there is one block with rank greater than 2. 

The next step, in the latter case, is to remove a second copy of k, corresponding to 
another A^-pair, and do this for all possible combinations of i and j. If for all 
this vanishes for arbitrary null k, then either r < m — 2 or r = m — 1 and there is one 
block with rank greater than 3. In this case, one has to go to the next step by removing 
one further copy of k. If on the contrary, for some pair there is a null k such that 
this does not vanish, then there must exist another I such that the same happens for 
the pair (i,l). Two possibilities arise here. If also the removal of the pair k^k 11 ' in (jSJ) 
gives a non-vanishing expression, then A A11 ... |tm is antisymmetric on the indices fa, fa, fa 
exclusively. Thus, r = m — 1 and A is a (m — l)-fold form with m — 3 one-index blocks 
and one block with rank 3. If on the other hand the removal of the pair k^k^ 1 in (j3J) 
vanishes for all null k, then there must exist a companion p for % such that the removal of 
the pair k^k^ p in (jSJ) is a vanishing expression too. Then A A11 ... A(m has antisymmetries on 
exactly the pairs fa, fa and fa;, fa, and no more antisymmetries. In this case r = m — 2 
and there are m — 4 one- index blocks and 2 blocks with rank 2. Observe also that, in 
general, the position of the antisymmetric indices is clearly identified. 

One has to proceed in this manner until the process stops or all the k pi for i — 1, . . . , m 
have been removed. In general, the principal directions are characterized by the particular 
null k's which annihilate the first non-generically vanishing expression of type (J2J) with 
several /c-suppressions. As is known, these principal directions can be characterized as 
the null k such that its inner product followed by the exterior product on all blocks of A 
vanishes, see 

(viii) Perhaps one could add to the previous list of properties the manifest and desir- 
able fact that the universal energy operator, or the derived mathematical energy tensor, is 
very easily remembered: just use one as defined in (|T} for each pair of corresponding 

indices of the tensor A® A; see ©. 

Let us consider some illustrative and important examples. 

Rank-2 tensors Let A^ u be an arbitrary rank-2 tensor. Then its mathematical energy 
tensor (J3J) reads 

Zapx^iA} = A aX Ap^ + Ap X A afJ , - g a fiA pX A p il - g XfJL A ap Ap p + ^g a pgx fl A pa A pa 

This can be written independently of any (anti)-symmetry properties of A. Suppose now 
that A pv happens to be a 2- form A pv = A^a = F^, such as the electromagnetic field. 
Its well-known energy-momentum tensor (which in this case is also the basic superenergy 
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tensor of F pu ) is 

T X ,{F} = F Xp F/ - \gx»F pa F»° 
so that the mathematical energy tensor of F pv can be rewritten now simply as 

£ap\fA,{F} = — 2F Q ,( A F A1 ) /3 — g aj3 T Xp {F] — gxfj,T a/3 {F} . 

This tensor satisfies the dominant property, as does T Xp {F}. Formula (J7J) specializes to 
£{a(3X[j,){F} = — 2g^ a /3Txf 1 ){F}. The energy density of F^ v can be equally computed with 
any of the two tensors due to 

S^ Xp {F}u a u^u x u p = 2T Xp {F}u x u^ 

for any unit timelike vector u. The principal null directions of F^ Vl defined by k p F p [ p k u ] = 
— that is, the null eigenvectors of F pv — , are equally well characterized by any of 

S a px„{F}^k x k^ = T x ,{F}k x k p = 0. 

Observe, however, that the mathematical energy tensor has more information than the 
energy-momentum tensor. For instance, the pairs of null vectors k, £ such that 

£ a ^{F}£ a fk x k p = 

define the null directions such that F pv £ p k v = — that is, all timelike 2-planes orthogonal 
to F^ v — , while these cannot be retrieved with T Xp {F}. 

Recall that any rank-2 tensor can be uniquely decomposed into its symmetric and 
antisymmetric parts, by means of A pu = At^ + A^ IU ] = S pv + F pv . Then from the 
definition of (jUJ) one derives 

£{A} = £{S} + £{F} + 2S{S, F} 

which allows us to put in relation the mathematical energy tensor of A with those of its 
symmetric and antisymmetric parts. Note that this was impossible with the superenergy 
tensors, as they had different ranks for S and F. Observe also that S Q/ 3 Xpj {S, F} = 
—£xfj,ai3{S, F} so that the (super)-energy density of A is the sum of those of S and F. 

Rank-3 tensors Let Ap^ v be an arbitrary rank-3 tensor. Then its mathematical energy 
tensor Q reads 

£a/3X[iTv{A} = A aXT Ap^ v + A^ Xr A apiV + A a ^ T A/3 Xu + Ap pT A aXv 
~ga/3 (A pXr A p ^ u + A ppr A p Xu ) — g Xp (A apT Af3 P u + A^ pr A a p u ) — g TU (A aXp Ap p p + Ap Xp A ail p ) 

+gapgX l iA p arA pa u + g apgrvA pXa A p ^ + gx^gruA apa Ap p g a pg X pg TU A pa ^A p 
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Suppose now that Ap^ happens to be antisymmetric on (say) the last two indices, so that 
A is in fact a double (l,2)-form App V = Ap^ = Lp^ u ; this is the case of the covariant 
derivative of the electromagnetic field, V pF^ v . The basic superenergy tensor of L is 



Formula (0) specializes to £ a p{\pTu){L} = -2T a p^p{L}g TU )- 

Consider finally the case where App V is a 3-form, that is to say, App V = A[pp U ] = Hpp U . 
The superenergy tensor of H (also called its energy-momentum tensor in this case) is 



—Qrv (H a \ p Hpp p + Hp\ p H a p p ) + 2g a pg\pT ru {H} + 2g a pg TU T\p{H} + 2g\pg TU T a p{H} 
and © reduces now to £{ a p\p TU ){H} = Qg {a pg XfJ T Tu) {H}. 



Comments analogous to that of rank-2 tensors concerning superenergy density and 
principal null directions are in order for all cases here. I omit them for the sake of brevity. 

Double (2,2)-forms One can do the same as before for rank-4 tensors, but in order to 
keep within a reasonable length let me just consider the important case of double (2,2)- 
forms, which include the Riemann and Weyl tensors, as well as the tensor product of 
the electromagnetic field with itself. Let A a p X p = A^^ by any double (2,2)-form, not 
necessarily symmetric in the interchange of the antisymmetric pairs. Its mathematical 
energy tensor is 





hence 





^gapg^s (Ap a \ T A p p V + Ap a p T A p \y) + g\pg TV (A a ^ pcr Aps p + Ap^ piJ A a s p ) 
+g-y5gTuT a p\p{A} + g^gxpTap^A] + g a pg\pT lSTU {A} + T-ySXp{A} 
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where 



T a f3\ l _ l {A} — A ap \ a Ap p ^ + A apfl(T Af3 p \ a — 

~n9aj3A pT \ r7 A p ^ — —gx^A api y T Ap p + —g a fig\ ll Ap Tav A p 
Z Z o 

is the basic superenergy tensor of A, also called its Bel tensor jU [TUJ 123 EH] (The 
mathematical energy tensor of any rank-4 tensor can be easily read off from the above 
formulas.) Formula ((Zj) becomes now 

£(a/3~fS)(\nTi;){A} = 2g^T 7 s)(\p{A} g TU ) . 

The principal null directions of A a p\^, defined by ^^k^A^p^k^ = 23], can be 
characterized by either T a p\n{A\k a kt ] k x k^ = or, equivalently, by 

&af3~f8\iLTi>{A}k k k^k k^k = . 

As remarked before, however, the mathematical energy tensor has more information than 
the superenergy tensor. For instance, the quadruplets of null vectors k, £, n, m such that 

S a ^sx„ T u{A}e a e (3 n"'n s k x k p m T m u = 

provide a richer geometrical structure for double (2,2)-forms. 

To finalize, recall that the Riemann tensor Rapxp of the Lorentzian manifold can be 
decomposed into the sum of the Weyl conformal tensor C a p\p, which is trace free, and a 
second part E a p\p depending on the Ricci tensor exclusively, see e.g. GDI 123 • This 
decomposition implies immediately 

£{R} = £{C} + £{E} + 2£{C, E}. 

which is a generic decomposition implying that of the Bel tensor into the Bel-Robinson 
and the so-called 'pure matter' and 'matter-gravity coupling' super-energy tensors [TUll25j . 

All in all, it seems that the definition is well-posed and mathematically sound, 
recovers all good properties of previous super-energy and energy-momentum tensors, is 
generic, valid in arbitrary dimension and for any tensor fields irrespective of their (anti- 
symmetry properties, easy to use, provides more information concerning the seed tensor, 
is potentially fruitful because many properties can be proved by using the generic formula 
(0) and the properties of G a/ 3 Xfl in (Q) rather than the cumbersome explicit formulas for 
super-energy tensors, and last but not least, it is easily remembered. 
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